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Wc present numerical results for chains of SU(2) BPS monopolcs constructed from Nahm data. The long 
chain limit reveals an asymmetric behavior transverse to the periodic direction, with the asymmetry becom- 
ing more pronounced at shorter separations. This analysis is motivated by a search for semiclassical finite 
temperature instantons in the 3D SU(2) Georgi-Glashow model, but it appears that in the periodic limit 
the instanton chains either have logarithmically divergent action or wash themselves out. 
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In this note we investigate numerically the properties of long chains of SU{2) BPS monopoles. 
In addition to the inherent mathematical interest in a search for such solutions, our physical 



motivation is to seek caloron (finite temperature instanton) solutions in the 3D SU{2) Georgi- 
Glashow (GG) model. At T = 0, such 3D GG instantons are mathematically identical to the 't 
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^ I Hooft-Polyakov monopoles [l|, |^ of a 3 -|- I dimensional Yang-Mills-Higgs theory; and at T > 0, 

calorons of the 3D GG model would be periodic in the Euclidean time direction, with period 
P = ^, and so would be "periodic monopoles". Recall that in 4D Yang-Mills theory, explicit finite 
action caloron solutions are known, both with trivial 0| and nontrivial P, 0, Q| holonomy along 
the Euclidean time direction. But in the 3D GG model, no finite action caloron solutions are 
known. This lack of caloron solutions is somewhat surprising for the following reason. At T = 
the 3D GG model is confining due to nonperturbative instanton effects P|. More recently it has 
been shown analytically Q], and confirmed by a lattice analysis |lo| . that the 3D SU{2) GG model 
has a deconfining phase transition in the Z2 universality class at nonzero T, and that the critical 
behavior at this transition involves both instanton and charged scalar degrees of freedom. It is 
therefore somewhat surprising that it has not been possible to construct finite action semiclassical 
solutions corresponding to finite T instantons. 
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To establish notation, consider the SU{2) Georgi-Glashow (GG) model in 3 dimensional space- 
time. This model is an SU{2) Yang-Mills gauge theory minimally coupled to a scalar field, <J>, in 
the adjoint representation, and with a symmetry breaking quartic scalar potential. The Euclidean 
action is 

A 



S = Jd^x ^tr + ^I)^$'^Z)^$'^ + - (^'^^'^ 



(1) 



Here = iA«r^ F^, = d^A, - d,A^ - ig[A^,A,], $ = and D^<^ = d^^ - ig[A^,^], 

where t'^ {a = 1, 2, 3) denotes the 2x2 Pauli matrix generators of su{2). Note that in 3 dimensional 
spacetime the coupling g has dimensions of (mass)^/^, as do the fields $ and A^. 3D GG instantons 
("monopoles") are solutions to the classical Euclidean equations of motion: 



D F 



ig[D,<^,<^] 



(2) 



No explicit periodic solutions to the instanton ("monopole") equations @ are known. The most di- 
rect way to study multi-monopoles is by ansatz, either radial or axial. The axial ansatz can produce 
chains of alternating monopoles and antimonopoles but does not appear to support chains of 



monopoles. Instead, here we restrict ourselves to the BPS limit (the scalar self-coupling A ^ 0), 
where there exist other approaches to constructing monopoles. A hint of periodic monopoles ap- 
pears in 



13], where doubly periodic 4D Yang-Mills instantons were studied - as one period shrinks 
one is left with 3D BPS monopoles periodic in the other direction. But the truly periodic limit is not 
fully understood. Direct analyses [l4| have found periodic BPS monopoles, but with divergent 
energy. Also, Bielawski has studied the n-monople problem for widely separated monopoles, 
showing that the moduli space is exponentially close to the hyperkahler Gibbons-Manton metric 

Q 

In the BPS limit the instanton/monopole equations (jSJ simplify to 
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In the BPS limit the energy/action reduces to 



d-^x tr K> 



'S'bps 



Thus, a convenient gauge invariant quantity to study is 



(3) 



(4) 
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In this paper we consider the Nahm construction of BPS monopoles, in which one can 

reconstruct the monopole fields ^ and from a solution of the associated Nahm equations, as 
described below. 




FIG. 1: An array of Dirac monopoles aligned along the X2 axis with separation /?. In the large J limit, 
within the central slab region — ^ < X2 < §, the magnetic potential goes over from the 3D l/r form to a 
2D hi{p) form within the region of radius of order J(3. 



It is useful to consider what one might expect to find. Outside a core region (whose scale we 
take to be f in this paper) a single BPS monopole behaves like an abelian Dirac monopole. Thus, 
one might expect (at least when the period is greater than the core size, which corresponds to low 
T in the 3D GG language) that a periodic BPS monopole would look like a periodic array of Dirac 
monopoles, outside of some core region around each monopole (se e Figure I). But if this were the 
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This is because for an array of 



case, the action/energy would be logarithmically divergent 
N = 2J + 1 Dirac monopoles, along the X2 axis with separation P, the magnetic field is B = — Vc/f), 
where the magnetic scalar potential is 

+J ^ 

n=-j ^xl + xl + {X2 - n(5Y 
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For any finite N, at truly asymptotic distances |x| ^ J [3, the scalar potential (p behaves as yj ^ y^, 
and the solution looks like a charge N Dirac monopole. But at intermediate distances, inside a 
sphere of radius of order J (3 (see Figure ^ the behavior is quite different. For a periodic array in 
the J — > oo limit we can restrict our attention to the central slab with — f < X2 < f , and write 

col I ^ ^\ 

(f) = constant + — j^^^^^=^^^^^= (7) 

n^oo \ ^xl + xl + {x2-nPy ^f^J 

Here the sum is extended to infinity and the divergence at the origin has been subtracted (the prime 
indicates that the n = term is excluded in the sum). Within this slab, at transverse distances 
inside some region of order Jf], (f) behaves as : ~ In (^ 2j+i ) ' "^bere p is the radial transverse 



distance: p = \Jx\ + x^. This illustrates the physical origin of the transition from the 3D Laplace 
Green's function to the logarithmic 2D Laplace Green's function. So, assuming the periodic BPS 
monopole is like an array of Dirac monopoles outside a core region, we might expect that the Higgs 
field magnitude |$| would behave as |$J_~ In p at large p. This leads to logarithmically divergent 



action/energy, as can be seen from 
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Is it possible to avoid this divergence? One possible solution is that the expectation of inde- 
pendent core regions is not realized. Since the BPS action/energy is proportional to the instan- 
ton/monopole charge, it is natural to expect that within the central slab region — ^ < 2:2 < § the 
net charge should be 1 even in the J — > 00 limit, and so it might be possible to define a large 
J limit in which the action/energy stays finite. This could happen if the core region spreads out 
significantly in the transverse directions. Our strategy here is simply to investigate numerically 
the large J limit of a class of monopoles which are manifestly BPS solutions for any J and which 
correspond to = 2 J + 1 monopoles equally spaced along the X2 axis. 

We use the Nahm construction of monopoles 3, ; which is based on first solving the " Nahm 
equations" , a set of coupled ordinary differential equations [rather than the coupled partial differ- 
ential equations in : 

dTi 1 

Here Ti(s), for i = 1, 2, 3, and s € [—1, 1], are N x N square matrices for an A^-monopole, satisfying 
three conditions: (i) Ti{s) is regular for s € (—1, 1), and has simple poles at s = ±1, with matrix 
residues defining irreducible representations of su{2); (ii) = — T^; and (iii) Tj(— s) = T^{s) (in 
some basis). It is not trivial to find solutions to © satisfying these conditions. Nevertheless, an 
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explicit solution is known which corresponds to a chain of N BPS monopoles equally spaced along 
a line 3, 3 



(this generalizes an explicit 2-monopole solution in 
n{s) = -yds(^(s + l)K A;^ i?i 
Ms) = -^nsf(s + l)K A;) R2 



Tsis) = -yCs((s + l)K 



k] R. 



(9) 



Here ds, ns and cs are elliptic functions, with real elliptic parameter < < 1, and K(A;) is 
the elliptic quarter period. The constant matrices Ri form an x irreducible representa- 
tion of su{2). For definiteness, we choose R3 = —2i diag ( J, J — 1, ... , — J + 1, — J), and Ri^2 to 



be antihermitean matrices with nonzero entries: {Ri)n,n+i = iRi)n+i,n = i\/n{2J — n + 1) and 
{R2)n,n+i = —{R2)n+i,n = —\/n{2J — 71+1). Since we are primarily interested in the central 
monopole in a long chain, we take = 2 J + 1 to be odd. 

Given such Nahm data, to construct the monopole fields ^ and which satisfy the BPS 
equation ©, one first solves the ID Weyl equation 



^l2N-^ - In ® X ■ T + iTj{s) ^ Tj^ V = , 



(10) 



where x S parameterizes the solution. There are two independent solutions, Va{s]x), [for 
a = 1, 2], of (fTfl|) that can be normalized as 



ds vl 



(11) 



Then the su{2) monopole fields are given in terms of the Va{s;x) by the following construction: 



^ab{x) = i J ^ds s vl{s; x) Vb{s; x) 



(12) 
(13) 



Af{x) = J ^dsvi{s;x)—Vb{s;x) 

An efficient numerical scheme for solving (|10|) has been presented in where symmetric A^ = 3 
and A^ = 4 monopoles were studied from their known Nahm data. We have used a similar numerical 
scheme, programmed in matlab. The main complication is that the Nahm data, which appears in 
(fTn|l . has singularities at the end points. Indeed, near s = +1 there are simple poles: 

i {Rl Ti + i?2 «) r2 + i?3 ® T3) 



iTj (s) (g) Tj ~ 



-iTj{s) (g) Tj ~ 



2 1 + s 

1 {Rl (g>Ti+ R2(g>T2 - Rs^ Ts) 

2 1 - s 



-1 
+1 



(14) 
(15) 
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These residues have eigenvalues J and — (J + 1), with degeneracies 2J + 2 and 2 J, respectively. The 
negative eigenvalues correspond to singular behavior of the vector v{s) at the end-points and so 
must be suppressed in order to satisfy the normalization conditions (lllj) . That is, as s — > =Fl) v must 
lie in the 2 J + 2 dimensional subspace spanned by the eigenvectors of the associated residue matrix 
with positive eigenvalues. These subspaces are independent of x and of the elliptic parameter k 
(which determines t he p eriod and so may be constructed once and for all for a given monopole 
number. Following 1^, we solve ()1U|) by shooting from either end of the interval, and matching 
at s = 0. The intersection of these shoots leaves just two independent solutions, vi and U2, each 
of which is a 2A^-component column vector, depending parametrically on x. Given vi and U2, we 
construct the magnitude squared of the Higgs field using (|12() . We focus on |<I>p since it is gauge 
invariant, its zeros indicate the locations of the monopoles, and because it determines the BPS 
action/energy density in (jl}. We implemented this numerical algorithm for three different values 
of the elliptic parameter: k = 0.99, k = 0.5 and k = 0.2, which we found correspond to period 
j3 = 3.25, (3 = 0.8 and (3 = 0.3, respectively, measured in units where the single-monopole core size 
is 1. 

First, consider the behavior of |$p along the direction of the monopole chain, which for the 
choice of Nahm data in (jUj) is the X2 direction. (A simple relabeling rotates the monopole in x 
space). Figure 1 shows along the X2 axis, for various values of the monopole number, and for 
k = 0.99, k = 0.5 and k = 0.2, respectively. For the k = 0.99 case one can clearly see the array of 
equally spaced zeros along the X2 axis. Note that the peak values of between adjacent zeros 
is very small, and also decreases as N increases. As N gets large vanishes within a region of 
length N(3^ outside of which it reverts to its A^-monopole form (compare with figure 1). But for 
k = 0.5 and k = 0.2 it looks as though |<l>p vanishes inside this region, —J/3 < X2 < J/3, already 
for J = 2. This is not the case, however, as can be seen from Figure 3, which plots the logarithm 
of |<I>p along the X2 axis. From these we can clearly see the periodic structure, with another two 
monopoles appearing on the X2 axis each time N is increased by 2. Notice how extremely small 
the peak values of become between adjacent monopoles as N increases. This is an indication 
that there are strong core interactions, preventing from getting anywhere near its asymptotic 
value between adjacent monopoles. We deduce the period to be /3 ~ 3.25 for k = 0.99, (3 ~ 0.8 
for k = 0.5, and [3 ~ 0.3 for k = 0.2. This period is measured in units of the single-monopole core 
size, so k = 0.99 corresponds to separated monopoles, which in the 3D Georgi-Glashow language 
corresponds to the low temperature limit. For k = 0.5 the separation and core size are comparable, 
while for k = 0.2 the separation is less than (roughly one third) the core size. In the notation of 



I*!^ along axis for k=0.99 
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FIG. 2: Plots of |$p along the X2 axis for elliptic parameter k — 0.99 (top - Fig 2a), k = 0.5 (middle - 
Fig 2b) and k ~ 0.2 (bottom - Fig 2c). The different curves correspond to increasing values of N, starting 
with iV = 1 at the top. For k = 0.99 one can see equally spaced zeros appearing along the X2 axis, but for 
k — 0.5 and k = 0.2 these zeros cannot be seen on this scale. The periodicity can be seen in Figure 3 which 
plots log(|$p). 
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FIG. 3: Plots of log(|(f>p) along the X2 axis for elliptic parameters k = 0.99 (top - Fig 3a), k = 0.5 (middle 
- Fig 3b) and k = 0.2 (bottom - Fig 3c). The minima clearly indicate the periodic structure along the X2 
axis, and one can count the increasing number of monopoles in the chain for successive curves, beginning 
with iV 1 at the top in each plot. 
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FIG. 4: Plots of |$p along the xi axis for elliptic parameter k = 0.99 (top - Fig 4a), k = 0.5 (middle - Fig 
4b) and fc = 0.2 (bottom - Fig 4c). The different curves correspond to increasing values of A^, starting with 
= 1 at the top. Note that for k ~ 0.5 and k ~ 0.2 these are very similar to the corresponding curves in 
Figure 2, but for k — 0.99 there is a marked difference from the corresponding curve in Figure 2. 
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FIG. 5: Plots of log(|(f>p) along the xi axis for elliptic parameters k = 0.99 (top - Fig 5a), k = 0.5 (middle 
- Fig 5b) and k ~ 0.2 (bottom - Fig 5c). Unlike the corresponding plots in Figure 3 along the X2 axis, there 
is no sign of any periodic structure in the xi direction. 
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FIG. 6: Plots of |$p along the 2:3 axis for elliptic parameter k — 0.99 (top - Fig 6a), k 0.5 (middle 
Fig 6b) and k ~ 0.2 (bottom - Fig 6c). The different curves within each plot correspond to increasing 
starting with = 1 at the top. 
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Ward's recent paper the cases k = 0.99, k = 0.5 and k = 0.2 correspond to small, medium and 
large values of the parameter C defining the ratio of the single-monopole core size to the separation. 

The plots in Figure 2 should be contrasted with the behavior along the xi axis, shown in Figure 
4. For k = 0.5 and k = 0.2 the plots are very similar to those along the X2 axis, with |<l>p vanishing 
inside a region of size of order Np. But for k = 0.99 the xi axis behavior is very different from that 
along the X2 axis. This is the limit where the individual monopoles are separated by a distance 
greater than the core size of a single monopole. To be sure there is not some small scale structure 
hidden in these plots, in Figure 5 we show log(|<&p) along the xi axis. In contrast to the plots in 
Figure 3, there is clearly no small scale structure along the xi direction. This also serves as a check 
that the extremely fine structure along the X2 axis is genuine, rather than a numerical artifact. 

The dependence along the X3 axis is shown in Figure 6. One might expect that in the large N 
limit the monopole solution within the central X2 period would become radially symmetric in the 
transverse {xi,xs) plane. However, this is not necessarily the case. For k = 0.99 the dependence 
along the xi and X3 axes is similar, as can be seen by comparing the first plot in each of Figures 
4 and 6, although there is somewhat flatter behavior near the origin along the xi axis. But for 
k = 0.5 and k = 0.2 there is a distinct asymmetry in this transverse plane. The asymmetry is more 
pronounced as N increases, and is also more pronounced for smaller values of k. Smaller k means 
smaller period so there is greater overlap, which presumably explains the increased asymmetry. 
Indeed, for k = the solution becomes axially symmetric in the (xi, X2) plane (with our convention 
of axes) P, [2^. This can be seen in Figs. 2c and 4c, which plot for k = 0.2 along the X2 and 
xi axes, respectively, and which are indistinguishable from one another on this scale, while being 
very different from the X3 behaviour in Fig. 6c. This should be contrasted with the k = 0.99 and 
k = 0.5 cases shown in Figs. 2a, b and Figs 4a, b, where there is a distinct difference between the 
behaviour of along the X2 and xi directions. 

Although all these plots are along the axes, we found them to be representative of what happens 
within a cubic region of size ~ (iV/?)'^ centred on the origin. Thus, we make the following empirical 
observations: 

1. In the large N limit, |<I>p vanishes along the xi and X2 axes within some region of order 
{N (3)"^ centred on the origin. This also implies that part of the action/energy density, {dl + f?!) 
vanishes within this region. 

2. In the large N limit, |<I>p exhibits three different behaviors along the X3 axis. At short 
distance the behavior is quadratic in X3, with a coefficient that decreases with N like j;^j2- At 
very long distances \^\ ~ 1 — j-^, as expected outside the core of an A^-monopole. But in the 
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intermediate region it is not yet clear what the dependence is. It can be fit weh in this regime by 
a linear behavior or by a logarithm. 

3. For period [5 much greater than 1 (the single-monopole core size) there is some remnant of 
symmetry in the transverse (xi,X3) plane, but this symmetry is badly broken for /? < 1. 

Could these structures, which are clearly equally spaced in the X2 direction for any finite A^, 
survive as finite action/energy periodic structures in the infinite limit? These numerical obser- 
vations show that there is considerable spreading of the core region in the transverse (xi, 2:3) plane. 
Moreover, for (3 <1 this appears to be almost one-dimensional. In this case, it might be possible 
to obtain a finite action/energy in the large limit. However, this would be at the price of the 
solution "washing itself out", because as V^|<^p spreads further out along the 3:3 direction with 
increasing N ^ correspondingly its peak value decreases. Thus, while this solution might appear to 
have a finite action/energy in the large limit, the action/energy density delocalizes and shrinks 
away to nothing. 

We reached a limit of simple numerical analysis in obtaining the plots reported here. It would 
be interesting now to investigate these empirical hints analytically at large A^, possibly using a 
semiclassical large spin J approach, and also by focussing directly on the (3^1 and [3-^1 limits. 
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